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ABSTRACT 

We show that the observed K velocities and periodicities of AM CVn can be recon- 
ciled given a mass ratio q « 0.22 and a secondary star with a modest magnetic field of 
surface strength i? ^ 1 T. We see that the new mass ratio implies that the secondary 
is most likely semi-degenerate. The effect of the field on the accretion disc structure 
is examined. The theory of precessing discs and resonant orbits is generalised to en- 
compass higher order resonances than 3:2 and shown to retain consistency with the 
new mass ratio. 
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1 INTRODUCTION 

, AM CVn is the prototype of a helium-rich class of ultra- 
' short period Cataclysmic Variable (CV) binaries. The eleven 
member systems consist of a white dwarf primary accreting 
. material, through Roche lobe overflow, from a companion 
' star that is itself degenerate or semi-degenerate. Nelemans 
, et al. (2001a) recognised two formation scenarios for these 
' systems. The "white dwarf family" arise from detached dou- 
ble degenerate white dwarf binaries that evolved into contact 
through gravitational wave radiation (Faulkner, Flannery & 
Warner 1972). The "helium star family" arise from systems 
I where a low-mass helium-burning secondary is brought into 
contact, again through gravitational wave radiation. Once 
mass transfer has reduced M2 < 0.2Mq, core helium burn- 
ing ceases and the star becomes semi-degenerate (Savonije, 
de Kool & van den Heuvel 1986; Iben & Tutukov 1991). 

Seven of the known AM CVn stars, including the proto- 
type, show periodicities in addition to their orbital modula- 
tions that are generally interpreted as arising from a precess- 
ing, non-axisymmetric, accretion disc. Several of the mem- 
bers show regular dips in brightness on a timescale of a 
few days. This suggests a similar phenomenon to the disc 
thermal instabilities in dwarf novae (Smak 1983; Cannizzo, 
Shafter k Wheeler 1988; Tsugawa & Osaki 1997) albeit with 
a helium rather than hydrogen dominated disc. Dwarf no- 
vae show regular outbursts where the luminosity increases 
by 2-5 magnitudes and a subset also show "superhumps" 
caused by a precessing accretion disc. Here, though, the de- 
fault state appears to be "high": ie. one in which helium 
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is ionized and the disc has a high viscosity. At some point 
the disc drops below a critical temperature, helium recom- 
bines and the disc switches to a low viscosity state. Mate- 
rial collects in the accretion disc and switches back to the 
high state at a second, critical temperature. These critical 
temperatures are normally converted to equivalent critical 
surface densities Ecrit. Exceptions are AM CVn itself, that 
has never been observed to dip, and CP Com which appears 
to be in a permanent low state (Warner 1995) . 

Nelemans, Steeghs & Groot (2001b) examined obser- 
vations of AM CVn in detail and arrived at a definite 
identification of Porb ~ 1028.73 s. Using the beat period, 
Pb = 13.38 h, standard precessing disc theory gives a mass 
ratio for the system q = 0.087. However, their measurements 
of the Hell 4286 A line suggested Ki = 53±6 km s~^ which, 
coupled with K2 = 210-280 km s~^, indicate a mass ratio 
in the range 0.19 < g < 0.25. The authors noted the incon- 
sistency implicit here and chose to rely on the commonly 
used mass-ratio inferred from disc theory. They also noted 
the alternative that the usual precession Pb-Porb-5 relation 
(given later in equation (12)) might be inappropriate. 

The purpose of this paper is to reconcile these two meth- 
ods of determining q by proposing that a modest secondary 
magnetic field is present in AM CVn. We suggest that the 
true AM CVn mass ratio is that derived directly from the 
K velocities and that the widely accepted Pb-Porb-9 rela- 
tion used by Nelemans, Steeghs & Groot (2001b) is indeed 
inappropriate in this case. 

Pearson, Wynn & King (1997) studied the effect of a 
secondary's field on the behaviour of dwarf novae: placing 
upper limits of Pli < 0. 16 T (/X2 < 7 x 10^* G cm^ ) in U Gem 
and Pli < 0.023 T (^2 £4 x 10^^ G cm^) in Z Cha. They 
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found that stronger fields had the effect of shrinking the disc 
radius and of exciting resonances in the disc. However, the 
disc was too small to access the 3:2 resonance normally ob- 
served in superhumping systems and instead higher order 
resonances appeared. Those higher resonances arc available, 
in principle, to all CV accretion discs that do not reach the 
required 3:2 resonant radius, even in the absence of a sec- 
ondary field. However, in practice, the tidal effect that drives 
the superhump resonance is unable to excite the higher or- 
ders before they arc damped out by viscosity. It is the ad- 
dition of the secondary field that enables these higher res- 
onances to be excited and their appearance or otherwise 
that enables us to constrain the secondary's magnetic field 
strength. We would expect a secondary with a magnetic field 
in AM CVn to have a similar effect: io. to shrink the disc, 
excite higher order resonances and to thereby alter the ap- 
plicable Pb-Povh-q relation. 

It has long been problematic to understand the origin 
of the magnetic field in the secondaries of CVs: particu- 
larly those systems that lie below the 2-3 hr "period gap" 
where the stars arc fully convcctivc. However, wo know that 
such a field must be possible, if not ubiquitous, from the 
fact that the Polars are able to lock in synchronous rotar 
tion through the interaction of the two stellar dipoles. This 
locking requires that the secondary has a magnetic moment 
fi2 ~ 10-^-^-fO-''^ G cm-'^ (King, Whitclmrst & Frank 1990) 
that is similar to that of the primary. Observationally, re- 
sults from a variety of open clusters (Jones, Fischer & Stauf- 
fer 1996; Stauffer et al. 1997; Torndrup ot al. 2000; Reid & 
Mahoncy 2000) suggest that the rate of magnetic braking 
does not reduce abruptly when stars become fully convec- 
tive {M <0.3Mq); implying the continued presence of a sig- 
nificant magnetic field. Such a field also has implications 
for the standard model of the period gap that relies on the 
rate of magnetic braking dropping sharply at the upper end 
(Andronov, Pinsonneault & Sills 2003). 

Recently, models have begun to emerge for the gener- 
ation of magnetic fields by fully convectivc, rapidly rotat- 
ing, young stars (Kiiker & Riidiger 1999) that may survive 
as relics into the main sequence lifetime (Kitchatinov, Jar- 
dine & Collier Cameron 2001). These fields are non-dipolar 
and highly non-axisymmetry. Such geometries, if repeated in 
fully convective CV secondaries, may have interesting con- 
sequences for the equilibrium alignment of Polars and the 
interpretation of the "hot spot" to which material is chan- 
nelled at the white dwarf surface. These models, however, 
have not so far extended into the extremely rapid rotation 
regime required for CV secondaries; especially with as short 
a Porb as AM CVn. We are also confronted, in this partic- 
ular case, with a secondary that is at least semi-degenerate 
with the unknown consequences for the field that such a 
structure and previous evolution will have. However, since 
for a dipolar field f oa oc r~^, it would require a rela- 
tively much weaker dipole moment to achieve the same local 
field strength in these short period, small separation systems 
than in wider long period CVs. It is, therefore, not unrea- 
sonable to expect to see effects in these systems that might 
not be apparent in their longer period counterparts. 

In section 2 we outline the important background the- 
ory to the problem; generalising the theory of superhump 
resonances to higher orders in section 2.1 and deriving com- 



ponent masses for the system in section 2.2. We describe 
simulations of the system in section 3. 



2 THEORETICAL BACKGROUND 

2.1 Predicting Superhump Periods 

Hirose & Osaki (1990) derived the general expression for the 
ratio of the disc precession Wp and orbital Worb frequencies 
in terms of the mass ratio and radius of disc material. Their 

equation (8) is 



Worb (l + q)^ 
where 

Bo{r) = ife'l 
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(Brumberg 1995) is the zeroth order Laplace coefficient 
given in terms of the hypergeometric function F, q = 
M2/M1 (< 1) is the mass-ratio and r is the radius of or- 
biting material expressed as a fraction of the separation a. 
Writing the hypergeometric function in its series form 
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The coefficients are calculable exactly but rapidly become 

unwieldy with large numerators and denominators. How- 
ever, after some algebra we can rewrite the general expres- 
sion for them as 



(l)(3)(5)...(2n-l) 

n!2" 

^(^)^ 



(5) 
(6) 



which at least clears away the less transparent fractional 
terms in Pochhammer's symbol of (|)^. Combining equa- 
tions (1), (4) and (6), and after some more algebra and a 
little differentiation, we can arrive at a general relation for 
resonant orbits 



Worb 4 

-I- c„r 

where the coefficients are related by 



(1 + g)^ 

,2(n-l) 
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(8) 



The values of C„ and Cn up to n = 5 are listed in Table 1 
for convenience. 

The radii of resonant orbits is given by Frank, King & 
Raine (1992). Restricting ourselves to the case of j:j — 1 
resonances, that are strongest, we obtain 



1 



i3 {l + qY 



(9) 
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n 





1 


2 


3 


4 


5 


Cn 


1 


1 

4 


9 

64 


25 
256 


1225 
16384 


3969 
65536 


Cn 





1 


16 


175 


3675 


72765 


8 


64 


1024 


16384 



Table 1. Coefficients contained in the expansions defined in equa- 
tions (6) and (8) up to n = 5. 



Resonance Beat Periods 




0.2 0.4 
Mass Ratio (q) 

Figure 1. Relation between the precession to orbital beat period 
and mass ratio given in equation (11) for different resonances 

j:j — 1. Also indicated is the range of the mass ratio allowed by 
K velocity measurements for AM CVn. 

Thus we finally arrive at 

LJr, 3 ( 



Worb 
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Since these are resonant orbits, the (long) beat period is 
related to the orbital and precession angular frequencies by 



Porb- 



OJorb 



J - 1 Wp 

This is plotted for several resonances in Figure 1. 
miliar approximation 



(11) 

The fa- 
(12) 



where A ^ 3.85 for 0.1 < g < 0.22 (Warner 1995), is recov- 
ered by setting j — 3 and evaluating the term in square 
brackets in (10) with q — 0.16. The limiting mass ratio 
q « 0.22 found by Whitehurst (1988a) arises from the largest 
value for which rs remains within the last stable stream line 
(Molnar & Kobulnickyl992). 

We can derive parallel approximations to (12) for other 
resonances in a similar way. The coefficient A need only be 
recalculated with the appropriate value of j. Table 2 gives 
A for different values of j; appropriate again to the range 
0.1 < q < 0.22 by evaluating with q = 0.16. 

For the particular case of AM CVn, using q = 0.22 the 
full expressions give rise to the predicted superhump periods 
and disc radii shown in Table 3. 



J 
A 



3 
3.86 



4 
5.33 



5 
6.75 



6 
8.16 



7 
9.55 



10.9 



9 
12.3 



Table 2. Coefficient A contained in the approximation (12) for 
several resonances and appropriate to 0.1 < g < 0.22. 



Resonance 
{3 -3-1) 


Radius 


^orb 


3 : 2 


0.450 


21.8 


4 : 3 


0.371 


29.8 


5 : 4 


0.320 


37.7 


6 : 5 


0.283 


45.5 


7 : 6 


0.256 


53.2 


8 : 7 


0.234 


60.8 


9 : N 


0.210 


()8. 1 



Table 3. The radius and beat period of superhump resonances 
for q = 0.22. 



2.2 Mass Determinations 

We can write the approximation by Eggleton (1983) for the 
size of the secondary's Roche lobe as 



R2 = 



47r2 



0.49g3(l -^g)3 
O.Gqi -|-ln(l-|-g! 



. !U"3 p3 
,-'"2 -^orb' 



(13) 



where we have eliminated a using Kepler's Third Law. 
Knowing Porb and q, for a given mass-radius relation we 
can eliminate R2 and solve for M2. 

For a fully degenerate He white dwarf secondary, the 
relation 



(10) ^ « 0.0106 - 0.0064 In f ^ ) -h 0.0015 f 



M 

Me 



M 
Me 



(14) 



(Zapolsky & Salpeter 1969; Rappaport & Joss 1984) is ap- 
propriate. For a semi-degenerate secondary, the relation is 
more problematic. However, approximations of the form 



R 
Re 



M 



(15) 



have been found by Tutukov & Fedorova (1989) with 6 = 
0.043, a = 0.062 (hereafter 'TF' parameters) and Savonije, 
de Kool & van den Heuvel (1986) with b = 0.029, a = 0.19 
(hereafter 'SKH'). The resulting values for Mi and M2 are 
listed in Table 4 for each case. 



Secondary Type 


Ml /Me 


M2/Me 


Degenerate 


0.156 


0.034 


Semi-degenerate (TF) 


0.529 


0.116 


Semi-degenerate (SKH) 


0.418 


0.092 



Table 4. Derived values for the AM CVn component masses 
assuming q = 0.22 and different mass-radius relations for the 
secondary. 
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Nelemans ct al. (2001a) carried out a population syn- 
thesis study ol AM CVn systems. They noted that the pre- 
viously adopted q = 0.087 gave rise to component masses 
that were difficult to reconcile with their results. For a de- 
generate secondary and that mass ratio, the system would 
lie at the low end of the predicted distribution and for a 
semi-degenerate system would imply a primary mass close 
to, if not in excess of, the Chandrasekhar mass. Locating 
the results from Table 4 on their Figure 1, we see that for 
q — 0.22 and a fully degenerate secondary, the system would 
have to have been born at an implausibly unlikely position 
in the tail of the distribution. For a semi-degenerate sec- 
ondary, however, the difficulties of a very high Mi have now 
been resolved with a more comfortable Mi ~ O.SMq. 



3 SIMULATIONS 

Modelling was carried out using the same code HYDISC 
(Whitehurst 1987; Whitehurst 1988b) used originally to 
demonstrate the tidal origin of superhump resonances 
(Whitehurst 1988a). It has since been adapted to a variety of 
other CV (Wynn & King 1995; Wynn, King & Horno 1997; 
King & Wynn 1999) and young star (Pearson & King 1995; 
Ultchin, Regev & Wynn 2002) applications. The modifica- 
tions necessary to include a secondary magnetic field were 
discussed in Pearson, Wynn & King (1997). 

Simulations of the effect of the secondary field in CVs 
were recently reexamined by Murray et al. (2002) using SPH. 
These authors modelled the magnetic interaction through an 
acceleration 



a= -fc[v-v/]±. 



(16) 



where v and v/ are the velocities of the material and field 
lines respectively and _L indicates that only components per- 
pendicular to the field linos are considered. Such a form, 
originating from models that used a diamagnetic prescrip- 
tion for the interaction, can also be applied to more general 
forms. In that case, our ignorance of the details of the inter- 
action are contained in k and, in particular, in its behaviour 
as a function of position and magnetic field strength. We 
prefer here to retain the model of the interaction contained 
in Pearson, Wynn & King (1997). That is, that the disc is 
assumed to wind up the magnetic field through advection 
until - B, (Aly & Kuijpers 1990; Wang 1996). This 
enables us to maintain an explicit relationship between the 
force and secondary field strength. Specifically, this earlier 
study showed that the acceleration experienced by each par- 
ticle could be modelled as 



_ 2ttB'^RAR 



(17) 



where rup is the mass represented by each particle, N is the 
number of particles in a bin of width AT? at a distance R from 
the white dwarf, B is the local magnetic field strength and /3 
is a dimensionless parameter related to the exact geometry 
of the field distortion (set equal to unity hereafter) . Since we 
have a strong constraint in the superhump period, we can 
place limits on the possible secondary field strength. 

The secondary's magnetic field acts on the disc mate- 
rial in a similar way to the gravitational tidal interaction. 
However, the disc material is now coupled more strongly 



to the secondary by the field. As a result, the disc needs 
to find a new configuration where the outward advection 
of angular momentum can meet the more efficient transfer 
back to the secondary. The B^ oc dependence provides a 
sharp "spike" to the interaction that acts to promote the ap- 
pearance of resonance phenomena. There is a limiting field 
strength below which the disc can find a suitable new struc- 
ture. Above this, angular momentum extraction is too effi- 
cient, the disc is completely disrupted and material is rapidly 
accreted by the primary. 

We adopted the system orbital period of 1028.7 s con- 
firmed by Nelemans, Steeghs & Groot (2001b), the mass 
ratio q = 0.22 in the middle of the range implied by the K 
velocities and Mi — 0.5Mq. 

Each simulation was started with a burst of particles 
over an initial 2Porb before being reduced to a steady 90 
per orbit. The magnetic interaction was "switched on" after 
3-Porb oiKx: a disc had begun to form. Each of the simula- 
tions was allowed to run for at least 500Porb, allowing an 
equilibrium to be reached. 

We can estimate the effective value of the Shakura- 
Sunyaev parameter a from the ring of material that forms 
initially at the circularization radius Rc in a non-magnetic 
simulation. Ignoring any additional mass transfer of mate- 
rial, the surface density E of an initial ring of mass m will 
evolve as 



exp< - 



Ii 



(t) 



(18) 



(Frank, King & Raine 1992) where x = R/Rc, / is a mod- 
ified Bessel function and t = VlvtR^^ . Fitting this to 
the early time behaviour where the much lower "steady 

state" mass transfer has negligible effect gives a viscosity 
J/ « 6.1 X 10** vc? s~^. From a Gaussian fitted to the vertical 
disc profile we have a scale height w 1.2 x 10® m. Hence, 
using a sound speed appropriate to a fully ionized gas of 
solar abundance at ~ 10'' K, we derive a ^ 0.3. 

The particle mass used in the magnetic interaction 
(17) was determined from a nominal mass transfer rate of 
10^* kg s~^. However, as the particles have negligible mass 
gravitationally, this is the only place in which rup appears 
in the equations of motion. Thus the magnetic simulations 
may be characterised by the parameter j^. We use the Li 
point as the fiducial position for the surface strength of the 
secondary's field and assume the field strength to vary as for 
a dipole. 

In common with Pearson, Wynn & King (1997) and 
Murray et al. (2002), we see that the magnetic field tends 
to shrink and promote resonant behaviour in the disc 
(see eg. Figure 2) although the resonances are relatively 
weak compared to SU UMa superhumps. Simulations with 
Sli £ 0.5 T are unable to excite any resonances. Similarly, 
when Bli ~ 1 T the disc is completely disrupted. Surface 
density profiles are shown in Figure 3 for different mag- 
netic field strengths. The critical surface densities for "dwarf 
nova" type thermal instabilities in helium discs are taken 
from Tsugawa & Osaki (1997) using a hot state a = 0.3 and 
cold state a = 0.03. We see that all the "magnetic" discs 
have surface densities between the two critical values and, 
hence, that they can maintain either a "hot" or "cold" state. 
However, the effective value of a derived above implies that 
the simulations represent a "hot" disc configuration which 
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Bli (T) 0.6 0.65 0.7 0.75 0.8 

28 35 43 51 58 

"j ' 1 5 () 7 8 



Surfoce Density Profiles 




7.2 7.4 7.6 7.8 

eg R (m) 

Figure 3. Mean surface density profile for equilibrium simula- 
tions with Sli = 0,0.6,0.65,0.7,0.75,0.8 T (decreasing radius). 
Also shown (dashed line) axe the critical surface densities for tran- 
sition to the "hot" (upper line) and "cold" states. 



is consistent with the observation that AM CVn is in a per- 
mament "high" state. 

We can derive the precession period of the disc by con- 
sidering plots of the angular distribution of particles, such 
as those in Figure 4, plotted every Porb- We can also use 
Fourier transforms of the number of particles in a particular 
azimuthal bin, although, even running for 512 Porb there is 
limited resolution at long P,. These periods are summarised 
for each field strength in Table 5. The models bracket the ob- 
served Pb = 46.8Porb and, conservatively, place limits on the 

/ . \-i/2 
allowed field strength of 0.65 T < Pli ( igii^tl s-i ) < 
0.8 T. These are equivalent to a magnetic moment, with the 
assumed mass transfer rate, of (6.0-7.3) x 10^^ G cm^ which 
is at the low end of the range deduced for Polar secondaries. 
We do not quote formal errors on these limits since they are 
dominated by the systematic effects in our assumptions re- 
garding the geometry of the field distortion. Pearson, Wynn 
&: King (1997) estimated that in a worst case these values 



Table 5. Measured beat period and closest predicted resonance 
for different field strengths. 

could be in error by no more than a factor 3 and were prob- 
ably significantly better. 



4 SUMMARY 

We have shown that the observed K velocities in AM CVn 
can be reconciled with our understanding of superhump res- 
onance behaviour given a modest secondary magnetic field. 
The field acts to shrink the disc and excite higher order res- 
onances than those tidally excited in SU UMa systems. The 
sensitivity of the disc behaviour to the magnetic field places 
a tight constraint upon the possible field strength of the sec- 
ondary although our errors are dominated by assumptions 
inherent in the model. A value of q — 0.22 strongly sug- 
gests that the secondary in AM CVn is semi-degenerate and 
is more easily reconciled with population synthesis results 
than the previous q = 0.087. 

We intend to return to this system in a future paper ex- 
ploring more exhaustively the possible parameter space (Mi, 
q, B, V etc.) which AM CVn can occupy and still reconcile 
the disc precession and K velocity measurements. Further 
direct measurements to confirm K\ would be extremely use- 
ful. 
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